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Abstract 

The structure of the mirror nuclei ^Be and is studied in a microscopic a + 
a+n and a+a+p three-cluster model using a fully antisymmetrized 9-nucleon 
wave function. The two-nucleon interaction includes central and spin-orbit 
components together with the Coulomb potential. The ground state of ^Be 
is obtained accurately with the stochastic variational method, while several 
particle-unbound states of both ^Be and ^B are investigated with the complex 
scaling method. The calculation for ^Be supports the recent identification for 
the existence of two broad states around 6.5 MeV, and predicts the and 
I2 states at about 4.5 MeV and 8 MeV, respectively. The similarity of the 
calculated spectra of ^Be and ^B enables one to identify unknown spins and 
parities of the ^B states. Available data on electromagnetic moments and 
elastic electron scatterings are reproduced very well. The enhancement of 
the El transition of the first excited state in ^Be is well accounted for. The 
calculated density of ^Be is found to reproduce the reaction cross section on 
a Carbon target. The analysis of the beta decay of ^Li to ^Be clearly shows 
that the wave function of ^Be must contain a small component that cannot be 
described by the simple a + a + n model. This small component can be well 
accounted for by extending a configuration space to include the distortion of 
the a-particle to t + p and h + n partitions. 

PACS number(s): 27.20.-Fn, 23.20.-g, 23.40.-s, 21.60.Gx 
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I. INTRODUCTION 



There has been a growing interest in the study of neutron-rich nuclei since the advent of 
radioactive nuclear beams. It was found [|l| that some light nuclei near the neutron drip-line 
exhibit neutron-halo structure or have thick neutron-skin clouds. The halo structure, a new 
form of the nuclear matter, is characterized by a spatially extended low density distribution 
around the core part of normal density. It is interesting to know how a nucleus changes its 
structure with the increase of the number of neutrons and how the binding of the neutrons 
is attained in such a system. In the very light nuclei the mean field is not stable enough to 
generate the regular shell structure but, instead, the clustering of the nucleons, especially the 
a-clustering plays an important role in determining their structure. Because of this the light 
nuclei show individual features which have strong dependence on the number of nucleons. 
The Be isotopes are of special interest in this respect because they show some anomalous 
features which are not easily understood in a simple shell model. Because the *Be nucleus 
is known to be a typical cluster state of two a-particles, it is interesting to attempt at 
describing heavier Be isotopes in a unified framework of two a-particles and extra neutrons. 
Our basic question is: How well does this picture give us a consistent understanding of the 
Be isotopes? This question naturally leads us to the application of a multicluster model. A 
fully microscopic multicluster model utilizes an A-nucleon wave function, incorporating the 
Pauli principle exactly. It has various applications in the structure study for the halo nuclei 
and in the nuclear astrophysics [^Q. 

The spectrum of ^Be is poorly known. This is probably because all the levels but the 
ground state are above the a + a + n threshold. Recent experiments have, however, 
isolated the broad level at 6.76 MeV to two states, the | , 6.38 MeV state and the | , 
6.76 MeV state. 

A few theoretical studies on ^Be have already been done in various models. A projected 
Hartree-Fock calculation |^ was carried out to study the electromagnetic properties of ^Be. 
A shell- model calculations in a {Q+l)Tiuj basis gave a reasonable spectrum but predicted 
too small dipole transition strength for the first excited state. There are several calcula- 
tions using an a -|- a -|- three-cluster model. Earlier calculations [|l3,|ll|] emphasized the 
three-body aspect of ^Be to explain its low-lying spectrum. These treated the a-particle 
as a structureless particle and considered its compositeness by redefining the potential with 
the Pauli correction. Recently, this type of macroscopic approach has been extensively ap- 
plied to the study of ^Be and ^B nuclei |T2|, by including the (aa)A^-type arrangement in 
the calculation. On the other hand, some microscopic cluster-model calculations starting 
from 9-nucleon wave functions were accomplished in the resonating group method |T^ 



or 



in the generator coordinate method ||1^,|T5|,|T6[ . Our microscopic multicluster model has the 
advantage that the distortion of the constituent clusters, e.g., the a-particle, when needed, 
can be included in the calculation in a consistent way. An example indicating this necessity 
will be discussed later in case of the (3 decay of ^Li to ^Be. The macroscopic model has, 
however, a difficulty in taking the possibility of the cluster distortion into account. 

The calculation of Ref. [|13| considered the three channels of ^Be(0"'")-|-n, ^He( g.s.)+a, 
and ^Be(2"'")-|-n to describe the levels of ^Be. A molecular model was applied in the generator 



coordinate framework to study the structure of Be [jT4|,|T5|. The calculation of Ref. |T6 



included only ^Be-|-A^ channel, where the motion of the two a-particles in ^Be was described 
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in a restricted space. The two generator coordinate method calculations gave a reasonable 
agreement with experiment. There are, however, some noticeable disagreements between the 
theory and experiment. Both of the macroscopic and microscopic calculations done so far 
were limited either in taking the Pauli principle into account or in treating the three-body 
dynamics. Further improvement will be attainable by treating the three-body dynamics 
more completely. 

As the first of the series of studies on the Be isotopes we show in this paper the results 
of calculation for ^Be in a microscopic a + a + n model. At the same time we consider its 
mirror nucleus ^B in an a -|- a -|- p three-cluster model. One of the main objectives in this 
paper is to assess the validity of our basic assumption in case of ^Be. This is substantially 
important for the study on heavier Be isotopes. To this end we carry out an extensive three- 
cluster model calculation that has no limitations mentioned above, and investigate carefully 
some important properties of the low-lying states of ^Be, that is, the energy spectrum, the 
magnetic and quadrupole moments of the | ground state, and the electron scattering form 
factors. Of particular interest is the enhancement of the electric dipole transition from the 
first excited state to the ground state |T^,[T8|. This reduced transition probability is 



nearly as large as the well known one of ^^Be. The mechanism of the enhancement in these 
cases may be related to each other. Another interest is the (3 decay from ^Li to the low-lying 
states of ^Be [ffOf- We will show that this /3 decay is useful to reveal small components 



contained in the wave function of ^Be. 

The levels of ^B are all particle-unbound and only few of them have spin assignments 
0. There are discussions on the missing i"*" state from the viewpoint of the Coulomb 
displacement energy [^,^]. As a mirror nucleus of ^Be, ^B can be described in an a + 



a+p three-cluster model. A cluster model has a unique advantage that it can describe the 
asymptotic part of a wave function well and thereby predict the position and width of a 
resonance. This is a very important ingredient for a detailed structure study of both ^Be 
and ^B because their states are mostly unbound. 

In our approach the total wave function is given as an antisymmetrized product of the in- 
ternal states of the clusters and the function of the relative motion. The antisymmetrization 
of all the nucleons is exactly taken into account. Two types of cluster arrangements, {aa)N 
and {aN)a, are combined to include the different correlation between the clusters. The nu- 
cleon in the {aa)N arrangement corresponds to moving in a "molecular" orbit around the 
^Be=(aa) core. On the contrary, the {aN)a arrangement is suited to describe an "atomic" 
orbit of the nucleon around the a-particle. This analogy should not, however, be taken 
so literally particularly when the particles come closer, because the configurations of the 
two arrangements have considerable overlap. The function of the relative motion is approx- 
imated by a linear combination of nodeless harmonic-oscillator functions of different size 
parameters. Our experience [|22| , |23| shows that the approximation with such functions gives 
an accurate description up to large distances. To keep the dimension of the basis low, we 
apply the stochastic variational method (SVM) [pSj^^ , in which we set up the "impor- 
tant" basis states stepwise by using an admittance test. This procedure was successfully 



applied to study the exotic nuclei and also to few-body systems ||23|] . 

The plan of this paper is as follows. In section II we give a brief outline of our formalism. 
The microscopic three-cluster model is presented in subsection IIA. The scaling methods 
which we apply to determine the position and width of a resonance state are briefly explained 
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in subsection IIB. Section III contains the results of calculations. The input parameters are 
given in subsection IIIA. The relative importance of the arrangements and the angular mo- 
mentum channels are discussed in subsection IIIB. Energies, radii, magnetic and quadrupole 
moments, electron scattering form factors are compared with experiment in subsection IIIC. 
The density distributions and the spectroscopic amplitudes are discussed in subsection HID. 
The (3 decay of ''Li to the states of ^Be is discussed in subsection HIE. In the last section 
we summarize the most important conclusions. 



II. FORMALISM 

A. A microscopic three-cluster model 

To describe the system consisting of a + a + n for ^Be or of a + a + p for ^B, we 
build up a trial function which is a sum over two cluster arrangements /i, /ii = {aQ)N 
and ii2 = {aN)a, with N = n or p. Each arrangement is associated with a particular set of 
intercluster Jacobi coordinates Pi and The coordinates Pi and P2 in the arrangement /ii 
are chosen to stand for the relative coordinate of the a-particles and the nucleon coordinate 
measured from the center-of-mass coordinate of two a's, while, in the arrangement fi2, they 
represent the relative distance vector between the nucleon and a and the relative coordinate 
of another a from the center-of-mass coordinate of the nucleon and a. The arrangement 
yUi is suited to describe the component corresponding to the ^Be-|-iV decomposition at large 
distances, while the arrangement //2 corresponds to the ^Re+a decomposition. The total 
orbital angular momentum L is obtained by coupling the orbital angular momenta ii = £f 
belonging to the Jacobi coordinates , and then it is coupled with the total spin 5* = | 
to get the total angular momentum J. See Fig. 1(a). The intrinsic wave function of the 
a-particle is constructed from a harmonic-oscillator Slater determinant with a fixed size 
parameter by eliminating the center-of-mass motion. The wave function of the intercluster 
motion is approximated by a linear combination of nodeless harmonic-oscillator functions 
(or "Gaussians" ) of different size parameters: 



^£m{^,p) = Ge{iy)ex.p{-i/p'^)yi^{p), 



(1) 



with 



2 2^+7/2 

0F(2^+1)!! 



1/2 



3^£m(x) = X%^{iL). 



(2) 



The wave function with the angular momenta [S, {iii2)L]JM {S =|) in the arrangement 
II can be written as 



K 



JM 



(3) 



where z/^, is the /cth size parameter of the ith relative motion in the cluster arrangement /i, A 
is the intercluster antisymmetrizer normalized such that the normalization kernel approaches 
the unit operator in the limit of infinite cluster separation, ^sms is a product of the intrinsic 
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wave functions of the two a-particles and the nucleon's spin-isospin function and K stands 
for the set of the indices {^1,^2} of the size parameters. By using an integral transformation 
p3| , the antisymmetrized product in Eq. (3) can be rewritten as a hnear combination of 
Slater determinants of Gaussian wave-packet single-particle functions. The matrix elements 
between Slater determinants of these nonorthogonal single-particle states are easily evaluated 
and can be expressed in a closed analytical form. 

The variational trial function is a combination of different arrangements and intercluster 
angular momenta: 

^JM= {^r5,(£i^2)L]JM + ^S(^i£2)L]Jm}- (4) 

It is noted that our wave function is fully antisymmetrized, free from the spurious center-of- 
mass motion (actually the total center-of-mass motion is eliminated) and has a good total 
angular momentum and parity. Our calculation is the so-called "variation after projection" 
type. 

The partial waves in a given cluster arrangement form a complete set of states and the 
different Jacobi coordinate systems are, therefore, equivalent in principle. One might thus 
think that we only need to choose a particular arrangement, /ii or /i2, and to decompose the 
wave function into a complete set of partial waves in this arrangement, and that the inclusion 
of both the arrangements implied in Eq. (4) would be redundant. Our experience shows. 



however, that the convergence of energy in a fixed arrangement is rather slow. The reason 

[s,{e[e'^)L'].jM 



is that the components ^(5 (£i^2)i]JA/ "^[s {£' e')l']jr4 arrangements /ii and fi2 are 



rather different, especially, at large distances and that any component ^[5 (£j^2)i]JM '^^'^ only 
be represented by an infinite sum in terms of the arrangement /i2. Moreover, the inclusion 
of high partial waves in the calculation is quite expensive. Our favorite choice is, therefore, 
to (1) decompose the wave function into partial waves in a given arrangement, (2) truncate 
the higher partial waves, and (3) complete the wave function by the inclusion of low partial 
waves of different arrangements. 

The arrangements and the angular momenta combined with the size parameters in the 
expansion make the dimension of the basis large. These basis functions are, however, 
nonorthogonal to each other and not all of them are equally important. In a previous paper 



25| we tested different methods to select the parameters z/^, that span most adequately the 



state space, while the dimension of the basis is kept feasible. The most efficient procedure 



found is the stochastic selection p2| , p3| : We generate size parameter sets by a random choice 
from a region which is physically important. The parameter sets that satisfy an admittance 
condition are chosen to generate basis states. The calculation was repeated several times to 
check the convergence. The dimension for the ^Be ground state is around 90. 



B. The scaling method for resonances 

Except for the ground state of ^Be, all the states of ^Be and ^B are above the three-body 
threshold. The I"*", 1.68 MeV state of ^Be lies just 111 keV above the threshold, but has 
a width of 217±10 keV. The |", 2.43 MeV state has a narrow width. The widths of other 
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states of isospin | range from several hundreds keV to about 1 MeV. The states of have 
generally wider widths than the corresponding states of ^Be. 

Resonances are associated with complex eigenvalues of the time-independent Schrodinger 
equation. It is not trivial to calculate the energy and the width of a resonance state for a 
complex system. Several methods have been developed to obtain these complex eigenvalues 
using square integrable functions. The most well-known methods are the complex scaling 
P7| and the stabilization [^] methods. 

The complex scaling method uses the unitary transformation which dilates the internal 
coordinates of the system according to x — >■ xe*^, making the resonant wave functions square 
integrable. The eigenvalues that are associated with metastable resonance states appear as 
such complex eigenvalues that are independent of the scaling angle 9, when it is larger than 
a critical angle, and the eigenvalues that are associated with nonresonant continuum states 



appear as complex eigenvalues which are dependent on the scaling angle |23]. One can 
expand the eigenfunctions of the complex scaled Hamiltonian in terms of square integrable 
basis functions as we did for bound states. The variation of the energy functional with 
respect to the trial function, however, yields a stationary rather than a minimum principle. 
Therefore, the stochastic basis selection procedure cannot be applied here, but instead, we 
will work on a basis with fixed nonlinear parameters. 



The stabilization method utilizes the discrete states calculated in a box of large size. 
The stabilization method can be combined with the stochastic variational method. In this 
case we select the basis parameters from a confined interval. 

These methods have been widely applied for two- and three-body resonances in atomic 
physics. Recently, nuclear physicists have also began to use the complex scaling method as 
a useful tool to locate two- [^] and three-body |0 resonances of nuclear systems. 



Due to the complexity of the problem both methods require extreme numerical accuracy. 
To be on the safe side, we used these methods only when they are certainly able to give reli- 
able results. That is, we used the stabilization method for narrow resonances, and calculated 
only the resonance energy because the calculation of the width would require an excessively 
large computational burden. For these quasibound states the stabilized wave functions can 
directly be used to calculate the matrix element of a physical operator because they are real. 
To locate wider resonances we used the complex scaling method. In this case we calculated 
both the width and the position. 

We have found that the energy of the narrow | state can well be obtained by diag- 
onalizing the Hamiltonian in a sufficiently large basis of Eq. (3). The resonance energy 
remains rather stable against the change of the basis set within a reasonable range. The 
wave function obtained in this way is used to calculate the electromagnetic transition rates. 
It is very difficult to do better than this because enclosing the wave function in a box as 
required by the stabilization method is not trivial for the three-body system. 

To apply the complex scaling method to the present case, we define the transformation 
U{9) which acts on the function of the intercluster Jacobi coordinates, pi and 

U{e)f{plp^,) = e'^'/'f{p1e^',p1e^'). (5) 

The eigenvalue problem of the transformed Hamiltonian Hq = U{6)HU{6)^^ is solved for 
each 9 value. A resonance state corresponds to an square integrable solution of the trans- 
formed Hamiltonian and may be described as in Eq. (4). When the basis function of Eq. (3) 
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is employed, the operation U{9)~^ on the relative motion function is equivalent to multiply- 
ing the size parameters, v'^^ and z/^^^ t>y e~^*^. The energy Er and the width of a resonance 
are obtained as the real and imaginary parts of a complex eigenvalue. Eg = E^ — I^F/j, of 
Hq, which remains unchanged for arbitrary values of 9 within an appropriate range. 



III. RESULTS 
A. Input parameters 

The internal state of the a-particle was approximated by Os harmonic-oscillator Slater 
determinant wave function of a size parameter i>=muj/2h. The value of u was chosen to be 
0.26 fm~^ to reproduce the experimental charge radius of the a-particle. The results are 
insensitive to the choice of the size parameter within a reasonable limit. 

We used Minnesota nucleon-nucleon interaction which is a sum of central and spin- 



orbit potentials of Gaussian form. The Coulomb potential was included. The strength of the 
spin-orbit force was taken from the set IV of Reichstein and Tang, which gives a good fit to 
N+a phase shifts. The central part of the Minnesota potential contains an exchange-mixture 
parameter u. The potential with u=l corresponds to a Serber type mixture. Decreasing the 
value of u from unity implies increasing repulsion in odd partial waves, while keeping the 
strength of even partial waves unchanged. It was set to u = 0.94 in order to reproduce the 
ground state energy of ^Be. The value of m = 0.94 is very close to the value (0.95) which is 
needed to well describe the a + a scattering in the resonating group method Thus our 
choice should give a realistic interaction between the a-particles. The value of m = 0.94 is, 
however, slightly smaller than the value of 0.97 recommended for the description of + a 
scattering. By fixing the u and u parameters as described above, the model contains no free 
parameter. No change of the potential parameters was made between ^Be and ^B. 



B. Cluster arrangements and angular momentum channels 

In our model the total spin is uniquely given by 5'=| so that the total orbital angular 
momentum can take either L = J— |orL = J+ |. Let us show that both values of L are 
needed by taking an example of the magnetic moment of ^Be. Quite probably (and this will 
be confirmed later) the orbital motion of the protons gives a moderate contribution to the 
magnetic moment of ^Be and only the spin part needs to be considered to get a reasonable 
estimate of the magnetic moment. The magnetic moment is then approximated by (J=|, 
L = 1 and L = 2) 




where Qsin) = —3.826 is the spin ^f-factor of the neutron in units of nuclear magneton and 
is the amplitude of the total orbital angular momentum L in the ground state wave function. 
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If the ground state is purely of L=l, then the magnetic moment becomes —1.913 which 
is in disagreement with the observed value of /iexp = —1.1778 fi^. An L =2 component of 
about 20 % admixture is needed to reproduce the observed value. We will see later that the 
potential chosen gives just the needed admixture. It is instructive to note that the magnetic 
moment for pure L=l case is equal to the Schmidt value of the single P3/2 neutron. 

Table I lists a set of arrangements and angular momenta used in the present calculation. 
We did several pilot calculations to know the relative importance of the arrangements and 
the angular momentum channels. When all the nine sets of Table I are used for the | 
ground state, the energy from the a + a + n threshold is obtained as —1.431 MeV and the 
root mean square (rms) radius of point nucleon is 2.50 fm. Let us call this a full calculation. 
When we exclude three sets belonging to the arrangement fi2 = {<yN)a with £1 = or 2, 
both energy and radius hardly change from the result of the full calculation; the overlap of 
the approximate wave function with the full wave function is 0.9995. This result is physically 
acceptable because the p wave is of prime importance for the interaction between the neutron 
and the a-particle. If we further exclude three sets belonging to the /i2 arrangement with 
ii = 1, then the energy increases to —0.32 MeV and the radius increases to 2.57 fm. This 
suggests that the arrangement fii = {aa)N (*Be+n-type configuration) alone is imperfect 
to describe the ground state even though the s and d waves are taken into account for the 
motion of the two a-particles. This consideration leads us to the remark that the calculations 
of Refs. [0,0 using only the ^Be+A^ channel should be accepted with some reservations. 
On the other hand, if we exclude three sets belonging to the /ii arrangement, then the result 
is very close to the full calculation; the energy loss is merely 34 keV and the overlap of the 
wave functions is 0.9991. We can thus conclude that the ^He+a-type configuration with 
£1=1 constitutes a very good approximation to the ground state wave function. As is seen 
from Table I, the angular momentum in the 112 arrangement is restricted to £1 = 1 for other 
states. 

For resonance states, particularly for high spin resonances the inclusion of high partial 
waves becomes important to obtain stable resonance parameters in the complex scaling 
method. The complex eigenvalue of the rotated Hamiltonian Hg is obtained by using the 
basis function of Eq. (3). The size parameters of the basis function are not selected randomly 
but are chosen as = uqp^'^ {k = 1, ...,K). The values of z/q and K are varied for each 
resonance to get stable values for its energy and width. The adopted value of K is about 10 
in the present calculation. The basis dimension used to diagonalize the rotated Hamiltonian 
is times the number of the sets listed in Table I. Figure 2 displays an example of the 
complex scaled spectra of ^Be for J'" = | and | . One can see, besides the discretized 
points corresponding to the three-body continuum, those points which lie on straight lines 
starting from the positions of the resonances of the subsystems. 



C. Energy spectrum and electromagnetic properties 

The calculated spectra of ^Be and ^B are compared with experiment in Fig. 3. The 
theoretical level sequence in ^Be has a good correspondence with the observed spectrum. 
The second | resonance is obtained at 4.3 MeV excitation energy. The other calculations 



l^ , [T5|Jl6| also predict the 1^ state. Although no such state is cited in Ref. 0, the calculated 
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resonance may correspond to the state at 5.59 MeV mentioned in Ref. 0. We get two broad 
overlapping resonances with | and |^ at about 6.5 MeV. This agrees with the conclusion 

of the recent experiments |^,^. We could not find a resonance with | around 8 MeV 
excitation energy in accordance with Refs. although such a state is parenthetically 

quoted in Ref. 0. Instead of this a | resonance is obtained at 7.9 MeV, which agrees with 
the result of Refs. [jl4|,|l5|. The spectrum of is less known experimentally compared to 
that of ^Be. The calculated spectrum is similar to the one of ^Be. We can predict the energy 
and the width of several resonances in ^B with the same accuracy as the case of ^Be. For 
example, our calculation predicts a missing | state at 2.43 MeV, which is in agreement 
with the result of a recent ^Be(p, n) reaction |3^ that located the | state at 2.83 MeV. 
Although no definitive spin assignment is made to the state at 2.788 MeV excitation energy 
0, our calculation supports a |^ assignment rather than 

Table II lists the energies and the widths of the unbound states calculated by the complex 
scaling method. The energies of the | states of both ^Be and ^B are in good agreement 
with experiment. Their widths, though extremely narrow, are reasonably reproduced; the 
calculated width of ^Be is about 2 times larger than the observed value, while the width of 
^B is about a half of the experiment. The calculation reproduces the widths of other states 
within a factor of 2. Our result is in better agreement with experiment than the calculation 



of Ref. 16 



There has been considerable effort to determine the location of the | state from the 



point of view of a Thomas-Ehrman shift |3^. We applied the complex scaling method to 
find a resonance with J'^ = ^ by including the arrangements and the angular momentum 
channels listed in Table I. The present calculation could not identify such a stable complex 
eigenvalue that can be interpreted as a resonance. To estimate the El transition strength, 
we increase the value of u to make the state particle-bound. 

The electromagnetic moments and the rms radii of proton, neutron, and nucleon, assum- 
ing pointlike nucleons, are included in Table III. Bare operators are used in the calculation. 
The charge radius of ^Be with the effect of the proton's finite size becomes 2.54 fm and fits 
the experimental value of 2.519 ±0.012 fm [^. The rms radius of neutron is larger than that 
of proton by 0.2 fm. Both the magnetic and the quadrupole moments of ^Be are reproduced 
very well. As was stated in subsection IIIB, the contribution of the proton's orbital motion 
to the magnetic moment is rather small (0.28 fi^) and the contribution of the spin part, 
— 1.45 fiN, corresponds to 15.1 % admixture of the L =2 component. The Ml and E2 tran- 
sition probabilities of the | state to the ground state are also well reproduced. The strong 
El transition of the |^ state is in reasonable agreement with experiment. The El transition 
strength depends on the description of the tail part of the wave function. With m = 1.0 
the energy of the state changes to 593 keV below the threshold and the B{E1) value 
becomes 0.24 W.u. in good agreement with experiment. With u = 0.98 the energy goes up 
to 206 keV below the threshold and the exterior part of the wave function that does not 
contribute to the transition grows, thereby reducing the B{E1) value to 0.18 W.u. To our 
best knowledge, this is the first theoretical calculation which has been able to reproduce the 
El transition probability in a consistent way. Reference [|T8| argues that the experimental 
El strength is enhanced to 0.38±0.07 W.u. if the unbound nature of the state is taken into 
account. 
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Table III includes the results of other models. The /i and Q moments of the shell 
model were determined by using an effective interaction which was chosen to reproduce both 
energies and static moments of Op-shell nuclei 0. These values are rather close to those of 
Cohen-Kurath (8-16) POT calculation |]3^,|5|. A shell-model calculation of (0-1-1)^0; model 



space cannot account for the enhancement of the B{E2) transition; with the effective 
charge of 0.35e it gives about one third of the experimental value. The El transition 
probability of the lowest state to the ground state was predicted to be only 0.03 W.u. 
0. Another shell- model calculation in a similar basis reproduces reasonably the B{E2) 
value by using a large effective charge for neutron, but again gives a very small B{E1) 
value. Although the calculation of Ref. |jT6|] using only ^Be+n channel gives a reasonable 
agreement with experiment, we have already pointed out that the ^He-|-a type configuration 



leads to further improvement. A similar remark applies to the calculation of Ref. |T2[, which 
indicates that the charge radius and the quadrupole moment are considerably smaller than 
experiment. 

Further test of the wave function of the ^Be ground state is performed by the electron 



scattering data pH]. The longitudinal electron scattering form factor is calculated in a 



first-order plane wave Born approximation through 

\FL{q)\' = ^ Y.\{^jf\\Mr\imj.)\\ (7) 

where Ze is the charge of the nucleus and the reduced matrix element of the operator T^[q) 
is defined by 

{JM\T^\J'M') = ) ' {J'M'kn\JM){J\\T''\\J'). (8) 
v2t/ + 1 

The charge density multipole operator M|^'^^(g) which occurs in the form factor is given as 
a function of momentum transfer q from the charge density operator 

MZ\q) = f Jiiqr)YUi) E ^—^^i^i " - r)dr, (9) 

i=i ^ 

where is the nucleon coordinate and Rem is the total center-of-mass coordinate. Note 
that our wave function contains no center-of-mass motion. Figure 4 compares the calculated 
charge form factor with the experiment P^J57|J5^. The correction for the finite proton size 



is taken into account by multiplying the form factor with the proton's form factor used 
in Ref. [^. Both monopole (CO) and quadrupole (C2) terms contribute to the charge 



form factor. No effort has so far been made to separate those contributions experimentally. 
Polarized electrons and targets will be needed to do such experiments. The agreement 
between theory and experiment is good. This is perhaps not very surprising because the 
present model reproduces both charge radius and quadrupole moment accurately. It is clear 
that the quadrupole deformation of the charge density is important at higher values. The 
deformation of the proton and neutron density distributions will be discussed in the next 
subsection. 

The transverse electron scattering form factor gives information on the nuclear current 
density. It is calculated from the expression 
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The symmetry consideration on parity and time reversal shows that the elastic form factor 
receives no contribution of the transverse electric multipoles of the current density j(r). The 
transverse magnetic multipoles are defined by 

TZ\<1) = jMqr)YT,,{r).j{r)dr. (11) 

Here the vector spherical harmonics are defined with unit vector e as 

= mr)e]i (12) 

and the current density consists of the convection and magnetization currents: 

j(r) = E ^ 2^^' ^Pi^i^i - Rem - r) +5(ri - Rem - r)pi| 

+ V X (^-^ J2 ^^i^i^i - ^cra ' r)^^^ . (13) 

Here p, is the momentum of the nucleon in the center-of-mass system and /ij is the magnetic 
moment of the nucleon in units of nuclear magneton. Figure 5 compares the calculated 
transverse form factor with the data of Refs. |5D,5T|. Both Ml and M3 contributions are 
important to get a satisfactory reproduction of experiment. Shell- model calculations 
needed a quenching factor of about 0.7 for the transverse form factors, while no quenching 
is needed in our model. We can conclude that the ground-state wave function of the present 
model reproduces consistently all the electromagnetic properties of ^Be. 



D. Density distributions and spectroscopic amplitudes 

The proton and the neutron density distributions, defined by 

A 

p(r) = (^jj I J2 '^(r. - Rem - r)P, | ^jj) = po(r) + J] p,(r)r,o(?), (W) 

(where Pi projects out the protons or neutrons) are also determined. For the ground state 
of ^Be we have monopole and quadrupole(£ = 2) densities. The density distributions, po{r) 
and p2{f), are shown in Figs. 6(a) and 6(b). They are related to the rms radius and the 
quadrupole moment as below: 

{r') = ^ / poir^dr, (15) 



Z Jo 

g = (1^)1/2 r p^^ry^dr. (16) 



An analogous relation can be defined for the neutron case. The quadrupole moment becomes 
5.13 fm^ for the proton and 3.86 fm^ for the neutron. The fact that the neutron quadrupole 
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moment is smaller than the proton quadrupole moment is understood by noting that the 
single neutron cluster moves between the two a-particles for the most of time and thus 
makes the neutron density less deformed. 

The 2.43 MeV, | and 6.38 MeV, | states of ^Be together with the ground state 
approximately follow a J{J + 1) rule and can be considered to form a rotational band 
with K=^ [Hil- /,From the experimental quadrupole moment of the ground state, the 
intrinsic quadrupole moment Qq of the band is estimated as 26.5 fm^ by using the relation 
Q = ( j+f)'[2j+3) ^0 ~ '^^^ transition probability within the band is related, in 

the collective model, to the Qq value by 

B{E2- KJ, ^ KJ2) = -^e'Ql{J,K20\J2Kf, (17) 

which predicts 23.9 e^fm^ and 9.98 e^fm^ for the | ^ | and | — > | transitions, 
respectively. The corresponding experimental values are 27.1±2.0 e^fm^ and 7.0±3.0 e^fm^ 
0. Since the collective model prediction agrees reasonably well with experiment, it may 
be possible to extract the intrinsic deformation parameter /Sq by using the relation /Sq = 
\/f theory gives /5o=0.89, which is close to the empirical deformation parameters of 

neighbouring nuclei, e.g., jSo ~1.13 for ^''Be and jSo ~0.82 for ^''C, while the corresponding Qq 
values are 22.9 fm^ and 25.0 fm^, respectively |^. The deformation parameter j3 associated 
with the density of Eq. (14) is estimated by assuming that it can be approximated by 
Ps{r/ (1 — + PY2o{r)) from a spherical shape Ps{r). The extracted value of /5 is close to 
1/5 of the Po value as expected by the collective model. 

The monopole densities of the protons and the neutrons may be used to calculate the 
total reaction cross section at high energies. It is given, in the Glauber theory |Q, as 



an = J[l- exp{-2Imx(b)}]t^b, (18) 

where b is the impact parameter and the phase shift function, x, is related to the densities 
of the target and the projectile through the thickness function, T(s) = / p(s, z)dz, by 

ix{h) = - jj Tp{s)TT{t)V{h + s - t)dsdt. (19) 

Here F is the profile function for the NN scattering. The monopole densities of the proton 
and the neutron were used to construct the density of ^Be. The value of ^Be for a 
Carbon target at 800 MeV/nucleon is calculated to be 850 mb with the parametrization of 
the profile function used in Ref. [^. The interaction cross section measured by Tanihata et 
al. [||] is not exactly the same as but approximately equal to the reaction cross section. Their 
value is 806±9 mb, which is in a fair agreement with theory. The reaction cross section of 
^Be on a Cu target was measured by Saint-Laurent et al. |^ at about 45 MeV/nucleon. 
They extracted the reduced strong absorption radius, tq ~ 1.13 fm, for ^Be by fitting their 
measured cross sections to the formula by Kox et al [^. This formula predicts = 825 ±20 
mb for the ^Be+^^C system at relativistic energies as listed in Table III. We again confirm 
that our density is reliable enough to reproduce the experiment. 

We showed in the previous subsection that the enhancement of the El transition of the 
first excited state in ^Be is reproduced well. To understand this we note that the El operator 
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is recast to ^ey^(R2 — Rat), where Hz and Rat are the center-of-mass coordinates of 
the protons and the neutrons, respectively. The enhancement of the transition should be 
therefore related to the excitation of the corresponding motion in the excited state. In the 
a + a + n model the El excitation is caused by the valence neutron. Figure 7 compares 
the monopole density of the state obtained with u = 1.0 with that of the ground state. 
The proton density becomes smaller near the center but reaches up to larger distances, 
indicating the increase in the mean distance between the two a-particles. The neutron 
density shows a substantial decrease at 1~2 fm and a significant increase beyond 3 fm. The 
proton and neutron rms radii increased from 2.39 to 2.94 fm and from 2.58 to 5.59 fm, 
respectively. Though the increase of the proton size is moderate, that of the neutron size 
is dramatic. The picture emerging from this analysis is the following: The valence neutron 
in the ground state is mostly confined between the two a-particles but, in the excited I"*" 
state, moves around them in a spatially wider region. It is easy to understand that the large 
El transition strength has naturally come out from the structure change of the underlying 
states. 

Another interesting quantity that helps to reveal information on the wave function is the 
spectroscopic amplitude which, in the angular momentum projected basis, is defined by 

9ti,i,)Lir,R) = ^(^{[$5[>^..(pO>^^.(P^)]l],m'^(Pi -^)^(P2 -i?)}l^JM). (20) 

Figures 8(a), 8(b), and 8(c) display the spectroscopic amplitudes of the ground state of ^Be 
for some channels of the arrangement fii = {aa)n, letting r and R represent the distances 
of the two a-particles and of the neutron from their center-of-mass, respectively. Some 
remarkable features are that all the amplitudes have a peak at r ~ 3.2 fm and R ~ 2.3 fm 
and /^-independent nodes at r = 1 fm (for the s wave between a's) and r = 2 fm (for the 
s and d waves). These characteristics are understood by the microscopic a-a cluster- model 
analysis for ^Be. The appearance of the nodes is understood in relation to the existence of 



the Pauli-forbidden states The norm of the amplitude, which is called the spectroscopic 
factor, becomes 1.03, 0.77, and 0.32 corresponding to three channels shown in Fig. 8. We 
note that the norm is different from the so-called amount of clustering. The amplitudes 
corresponding to the arrangement fi2 = {an)a are plotted in Fig. 9, where r is now the 
distance between n and a and R the distance between their center-of-mass and a. The 
nodes appear also in this case but their positions alter particularly at large r. This is due to 
the fact that R is approximately equal to the a-a distance at small r but deviates largely 
from it with increasing r. The spectroscopic factor is 0.84 and 0.61, respectively. 



E. Beta decay of the ^Li ground state to ^Be 

Because the ground state and the | , 2.43 MeV state are described well by the present 
model, the P decay of the ^Li ground state to these states is expected to further test the 
accuracy of their wave functions or an available wave function of ^Li. The experimental 
value of log/t for the /3 decay to the ^Be ground state is about 5.31 |]7|JT9[|, indicating that 
the /3-decay matrix element is fairly suppressed despite the allowed transition. The weak 
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(5 decay is ascribed to the fact that the spatial symmetry of the main component of ^Be is 
different from that of ^Li |^|. The Gamow- Teller (GT) matrix element, 



M< 



f) 



9 
k=l 



(21) 



to any state of ^Be, if it is described by the a + a + n three-cluster model, always vanishes 
regardless of the wave function of ^Li. This is most easily understood by acting the Hermitian 
conjugate of the GT operator on the ^Be wave function and by noting that the spin-isospin 
part of the a-particle wave function is fully occupied. 

The above discussion indicates that the simple three-cluster model for ^Be must be 
modified to explain the /5 decay in spite of the successful results obtained in the previous 
subsections. The modification must be small enough not to destroy the agreement between 
experiment and the three-cluster model calculation. One possible way for the modification 
is to introduce the distortion of the a-particle into t+p and h + n configurations. To explore 
the consequence of this modification, let us assume that the intrinsic wave function of the 
a-particle can be expressed by 



(pa = VT 



+ CC 



(e) 



(22) 



where 0^°^ represents the part which can be described by the Os harmonic-oscillator Slater 
determinant, while (\)^^ the distorted part which is orthogonal to The ^Be wave function 
in a more realistic three-cluster model can therefore be approximated by 

Iv^j/^Be)) =M{{\- c2)|vl/W(9Be)) + 20^1^1 ^W^i^)^) + c^\4>f 4>t^ n)) . (23) 



Here the normalization constant. A/", is close to unity when c is small and it is suppressed 
below. The first term, |\l/j°''(^Be)) = \(\)^^ (\)^^ n) ^ is nothing but the one described by the 
a -|- a -|- n model and has no contribution to the /5 decay. By neglecting the last term, the 
GT matrix element is given by 



MGT{i ^ f) = 2cVl - 



n 



k=l 



2(v^(°H'Be)|vl/',^,(^Be))(vl/',^.(9Be) || Y.t-(k)a{k) ^ j^'U) 



k=l 



with 



l^'.^^Be)) 



(24) 



(25) 



Equations (24) and (25) are our basic equations to calculate the [3 decay matrix element 
when the distortion of the a-particle is included. 

The wave function of eq. (25) is obtained by extending the three-cluster model to the 
four-cluster model which includes a + t + p + n and a + h + n + n partitions. In order 
to avoid excessive numerical calculations, the angular momentum channels and the cluster 
arrangements are rather limited. See Fig. 1(b) and Table IV. The spins of all the clusters 
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were coupled to 5* = |. The isospins were not coupled to a definite value so that the wave 
function of the extended model may in general contain the total isospin of T = ^ and |. 
The potential favors T = ^ for the low-lying states of ^Be. 

The intrinsic wave function of the t- and /i-cluster was described by Os harmonic-oscillator 
Slater determinant of the same size parameter u as that of the a-particle. The ground state 
wave function obtained in the four-cluster model using the Minnesota potential of u=0.94 
has the overlap integral of 0.971 with the one obtained in the three-cluster model. Therefore, 
this new wave function should yield substantially the same results as the previous one for 
the electromagnetic properties. This is just what we have expected to maintain in extending 
the model space. 

To calculate the /3-decay probability we use the ^Li ground-state wave function which 
was obtained in a microscopic a + t + n + n four-cluster model |]22|. This model for ^Li 



reproduced both magnetic and quadrupole moments of the ground state very well. To fit 
the energy of the ^Be ground state to its experimental value from the four-body threshold, we 
changed the u parameter to 0.88 in the four-cluster model calculation. The overlap integral 
of the wave functions between the three-cluster and four-cluster models becomes 0.973 and 
the resulting log/t value is 5.60. The log/t value is still a little too large compared to the 
experimental value of 5.31, but this calculation strongly indicates that we are on the right 
track. A further refined four-cluster model calculation for both ^Be and ^Li will reduce 
the disagreement between experiment and theory because such a calculation is expected to 
enhance the GT matrix element. The shell-model calculation gives the log/t value in 
the range of 4.86 to 5.64, depending on the interaction used. It is interesting to analyze in 
the way presented above the /3 decay of to the low-lying states of ^B and an asymmetry 



in the /3-decay matrix elements of ^4 = 9 nuclei [49,19 



IV. SUMMARY 

The microscopic multicluster model was applied to the study of the mirror nuclei ^Be 
and ^B. They were described in a three-cluster model comprising two a-particles and a sin- 
gle nucleon. The two-nucleon interaction consists of the central and spin-orbit potentials 
together with the Coulomb potential. The same two-nucleon potential was employed for 
both ^Be and ^B. The ground state of ^Be, an only particle-bound state in this study, was 
obtained with the stochastic variational method, while the other particle-unbound states 
were studied by the complex scaling and the stabilization methods. The three-body dynam- 
ics of the clusters was taken into account by including both of the arrangements, {aa)N 
and {aN)a, and by using relevant partial waves between the relative motion of the clusters. 
The calculated spectrum of ^Be below an excitation energy of 8 MeV was in fair agreement 
with experiment. We obtained two broad overlapping resonances with J'^ = | and |^ 
around 6.5 MeV, in agreement with the conclusion of the recent experiments. Two states, 
I2 and I2 , were predicted at about 4.5 MeV and 8MeV in excitation energy, respectively. 
The spectrum of ^B was found to be similar to that of ^Be. The spin and parity of several 
states of ^B were predicted. The first excited I"*" state was not localized in the present study 
and thus no definite argument was possible on a Thomas-Ehrman shift in this case. 

The theory reproduced very well the electromagnetic properties of the ^Be ground state 
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such as the charge radius, the magnetic moment, the quadrupole moment, and the clastic 
electron scattering form factors. The calculated ground state density was consistent with 
the total reaction cross section data. The intrinsic deformation parameter of the density 
was found to be 0.89. The ^'^ |" El transition and the |~ |~ E2/M1 transitions 
were studied by treating the excited states as bound. The calculated transition rates were 
in good agreement with experiment. 

The fact that the present calculation reproduced all the data well strongly supports 
that the three-cluster model is quite appropriate for describing the structure of ^Be and 
^B, provided a full account of the dynamics is taken into account in the calculation. We 
were also able to understand the P decay of ^Li to ^Be by admixing the small components 
that are induced by the distortion of the a-particle into t + p and h + n configurations. A 
unique advantage of the microscopic multicluster model was exemplified by being able to 
accommodate such distortion in the model consistently. The study on heavier Be isotopes is 
in progress in the framework of the microscopic multicluster model including two a-particles 
and several neutrons. 

This work was supported by a Grant-in Aid for Scientific Research (No. 05243102 and 
No. 0664038) of the Ministry of Education, Science and Culture (Japan) and by OTKA 
Grant No. T17298 (Hungary). Most of the calculations were done with the use of RIKEN's 
VPP500 computer. 
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TABLES 



TABLE I. A set of arrangements and angular momenta included in the three-cluster model 
calculation for ^Be {N = n) and ^B {N = p). See Fig. 1(a) for the angular momenta £i and £2- 



J'^ arrangement angular momentum (£1, £2) -f' 



1/2- 


{aa)N 
(n'N)ry 


(0,1)1 
M Oil 


(2,1)1 

i^ 211 


(2,3)1 








1/2+ 


{aa)N 


(0,0)0 

n 1^0 


(2,2)0 

(1.1)1 


(2,2)1 








3/2- 


{aa)N 
laN)a 


(0,1)1 
(0,1)1 


(2,1)1 
(1,0)1 


(2,1)2 
(2,1)1 


(1,2)1 


(2,1)2 


(1,2)2 


3/2+ 

/ 


(aa)N 
{aN)a 


(2,2)1 
(1,1)1 


(0,2)2 
(1,1)2 


(2,0)2 
(1,3)2 


(2,2)2 


(2,4)2 


(4,2)2 


5/2- 


{aa)N 
laN)a 


(2.1) 2 

(1.2) 2 


(2,3)2 
(1,2)3 


(0,3)3 


(2,1)3 


(2,3)3 




5/2+ 


{aa)N 
laN)a 


(0,2)2 
(1,1)2 


(2,0)2 
(1,3)2 


(2.2) 2 

(1.3) 3 


(2,2)3 






7/2- 


{aa)N 
{aN)a 


(2.1) 3 

(1.2) 3 


(0,3)3 
(1,4)3 


(2.3) 3 

(1.4) 4 


(4,1)3 


(2,3)4 


(4,1)4 


9/2+ 


{aa)N 
{aN)a 


(2.2) 4 
(2,4)5 

(1.3) 4 


(0,4)4 
(4,2)5 
(1,5)4 


(4,0)4 

(4.4) 5 

(1.5) 5 


(2,4)4 


(4,2)4 


(4,4)4 
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TABLE II. Energies and widths of the unbound states in ^Be and ^B. The energy is from the 
three-body threshold. The spin and parity of the 3.065 MeV state of ^B is assumed to be I"*". 
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0.30 


0.004 


1/2+ 
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-0.7 






5/2- 


2.638±5 


0.081+5 


2.55 


0.044 


1/2- 


3.11" 


3.1 


2.73 


1.0 


5/2+ 


3.065±30 


0.550+40 


3.5 


1.2 


3/2+ 






4.6 


2.7 


3/22- 






4.2 


1.4 


7/2- 


7.25±60 


2.0+200 


7.0 


1.7 


9/2+ 






6.6 


3.3 


5/22- 






8.4 


2.4 



a) Ref. 0. b) Ref. 



c) Ref. §2 



20 



TABLE III. Radii and electromagnetic properties of ^Be. The reduced matrix elements are 
given in Weisskopf units. The bare-nucleon charges and (7-factors are used in the present calculation. 
The effective charges were used in the shell model calculation of Refs. and P to calculate the 
quadrupole moment and the E2 strength. See text for the B{E1) value of the present calculation. 







exp.'^ 


present 


Ref. [15| Ref. 16] Ref. 5| 


Ref. [|] 


3/2" 


E (MeV) 

Tm (fm) 

Tp (fm) 

Tn (fm) 


-1.5735 
2.37±0.01 


-1.431 
2.50 
2.39 
2.58 


2.62 


-0.89 








(/^jv) 
Q (e fm2) 


-1.1778±0.0009 


-1.169 


-1.23 


-1.52 


-1.27 


-1.070 




5.3±0.3 


5.13 


5.76 


4.77 


4.35 


4.66 




an (mb) 


825±20'' 


850 










5/2- 


E (MeV) 


0.8559 


0.883 




1.89 








5(^2; f^f) 


24.4±1.8 


22.0 


24.7 


23.5 


12.5 


~ 7 




E(M1;|--.|-; 


0.30±0.03 


0.229 




0.10 


0.23 





1/2+ E (MeV) 

B{E1; l^^f 



0.111 
0.22±0.09 



0.24 
0.18 



0.05 
0.68 



0.03 



0.75 
0.03 



a) Ref. [0. b) Ref. [|§ 
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TABLE IV. A set of arrangements and angular momenta included in the four-cluster model 
calculation for the '^Be ground state. Sec Fig. 1(b) for the angular momenta ii, £2, and £3. The 
spin of the nucleon clusters is coupled to S23. The total spin S is restricted to ^. 



J'^ arrangement 


angular momentum 




S23 


Q"r> /n— //j. \ \ 

^Be: 3/2 [[tp)a)n 


[(0,0)0,1]1 


[(0,2)2,1]1 


[(0,2)2, 1]2 


1 


{tp){an) 


[(0,0)0,1]1 


[(0,2)2,1]1 


[(0,2)2,1]2 


1 


{(tp)n)a 


[(0,1)1,0]1 


[(0,1)1,2]1 


[(0,1)1,2]2 


1 


{(tn)a)p 


[(0,0)0,1]1 


[(0,2)2,1]1 


[(0,2)2,1]2 


1 


{tn) {ap) 


[(0,0)0,1]1 


[(0,2)2,1]1 


[(0,2)2,1]2 


1 


{{tn)p)a 


[(0,1)1,0]1 


[(0,1)1,2]1 


[(0,1)1,2]2 


1 


{{hn)a)n 


[(0,0)0,1]1 


[(0,2)2,1]1 


[(0,2)2, 1]2 





{hn){an) 


[(0,0)0,1]1 


[(0,2)2,1]1 


[(0,2)2,1]2 





{{hn)n)a 


[(0,1)1,0]1 


[(0,1)1,2]1 


[(0,1)1,2]2 
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Figure Captions 



Fig. 1. Different arrangements used in the three-body (a) and four-body (b) calculations. The 
small circles are nucleons, the medium-size circle is a-particle, and the gray circle is 
3iV-cluster, t or h. The orbital angular momenta for the relative motion between the 
clusters connected by solid lines are denoted by £j. The spin of the clusters is Sj = |; 
the spin of the a-particle is zero and it is omitted. 

Fig. 2. Complex eigenvalues for J'^ = | (a) and | (b) of ^Be. The rotation angle 9 is in 
units of radian. The point indicated by an open circle corresponds to a resonance. 

Fig. 3. Experimental and calculated energies of ^Be (a) and ^B (b) from the three-body thresh- 
old. The data are from Refs. PJ^J^l- "^^^ 3.065 MeV state of ^B is assumed to be 



Fig. 4. Elastic charge form factor for Be. The data are from Refs. [ p6| , p^ ,p8[] . 



Fig. 5. Elastic transverse form factor for ^Be. The data are from Refs. POipTI . 



Fig. 6. Monopole (a) and quadrupole (b) density distributions of protons and neutrons for 
the ground state of ^Be. 

Fig. 7. Monopole density distributions of protons and neutrons, (a) in linear scale and (b) in 
logarithmic scale, for the excited | state and the | ground state of ^Be. The value 
of M = 1.0 is used for the I"*" state. 

Fig. 8. Spectroscopic amplitudes of the ground state of ^Be for the ^Be-|-n arrangement. The 
symbols r and R denote the distances of two a-particles and of n from their center-of- 
mass. The set of angular momenta is (a) £i = 0, £2 = 1, = 1, (b) £1 = 2, £2 = 1, = 1, 
and (c) i^ = 2j2 = l,L = 2. 

Fig. 9. Spectroscopic amplitudes of the ground state of ^Be for the ^He-t-a arrangement. The 
symbol r is the distance between n and a and R the distance between their center-of- 
mass and a. The set of angular momenta is (a) £1 = 1, £2 = 0, L = 1, (b) ii = 1, £2 = 
2, L = 1, and (c) £^ = 1^2 = 2,1 = 2. 
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(a) 3 -body cluster arrangement 
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(b) 4-body cluster arrangement 
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